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ABSTRACT 

By performing three-dimensional (3D) simulations that demonstrate the 
neutrino-driven core-collapse supernovae aided by the standing accretion shock 
instability (SASI), we study how the spiral modes of the SASI can have impacts 
on the properties of the gravitational-wave (GW) emission. To see the effects 
of rotation in the non-linear postbounce phase, we give a uniform rotation on 
the flow advecting from the outer boundary of the iron core, whose specific an- 
gular momentum is assumed to agree with recent stellar evolution models. We 
compute fifteen 3D models in which the initial angular momentum as well as the 
input neutrino luminosities from the protoneutron star are changed in a system- 
atic manner. By performing a ray-tracing analysis, we accurately estimate the 
GW amplitudes generated by anisotropic neutrino emission. Our results show 
that the gravitational waveforms from neutrinos in models that include rotation 
exhibit a common feature otherwise they vary much more stochastically in the 
absence of rotation. The breaking of the stochasticity stems from the excess of 
the neutrino emission parallel to the spin axis. This is because the compression 
of matter is more enhanced in the vicinity of the equatorial plane due to the 
growth of the spiral SASI modes, leading to the formation of the spiral flows 
circulating around the spin axis with higher temperatures. We point out that 
a recently proposed future space interferometers like Fabry-Perot type DECIGO 
would permit the detection of these signals for a Galactic supernova. 



Subject headings: supernovae: collapse — gravitational waves — neutrinos - 
hydrodynamics 
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Introduction 



The successful detectio n of neutrinos from SN1987A paved the way for Neutrino As- 
tronomy (IHirata et al.lll987l ). an alternative to conventional astronomy by electromagnetic 
waves. Core-collapse supernovae are now expected to be opening yet another astronomy, 
Gravitational- Wave Astronomy. Currently long-baseline l aser interferometers such as LIGO 
(Abbott et al.ll2005h . VIRGO 1 , GEO600 2 , and TAMA 300 (Undo fc the TAMA Collaboration 
20051 ) are operational (see, e.g., iHough et al.l (120051 ) for a recent review). For these detec- 
tors, core-collapse supernovae have been prop osed as one of the m o st pl a usible sources of 
gravitational waves (GWs) (see, for example, iKotake et al.l (120061 ); lOttl (120091 ) for recent 



reviews . 



Although the mechanism of explosion is not completely understood yet, current multi- 
dimensional simulations based on refined numerical models show several promising scenarios. 
Among the candidates is the neutrino heati ng mechanism ai d ed by convection and stand - 



ing accretion shoc k instability (SASI) (e.g., iMarek &: Janka (120091) ; iBruenn et al.l ( 120101 ); 



Suwa et al.l (120091 )). the acoustic mechanism dBurrows et al. 



2006h . or the magnet o hydro 



dynamic (MHD) mec hanis m fe.g.. iTakiwaki et al.l ( 12004 . 120091 ); lObergaulinger et al.l ( 120061 ); 
Burrows et al.l ( 120071 ). and IKotake et al.l (120061 ) for collective references therein). For the 
former two to be the case, the explosion geometry is expected to be unipolar and bipolar, 
and for the MHD mechanism to be bipolar. Since the GW signatures imprint a live infor- 
mation of the asphericity at the moment of explosion, they could provide us invaluable clues 
to understand the supernova mechanism. 

Traditi onally, most of the t heore ti cal predictions of GWs h ave focused on the bounc e 



signals (e.g..lMonchmever et al.l (ll99ll) ; IZwerger fc Mullerl (ll997h:lKotake et al.l ( 
Shibata fc Sekiguchil (boO^llott et all J20p4 12007 J ch: lDimmelmeier et al.l f|2002 



2003 



2007 



20041 : 



20081 ) 



Scheidegger et al.l (120101 )) and references therein). However recent stellar evolution calcula- 
tions suggest that rapid rotation assumed in most of the previous studies is not canonical 
for progenitors with neutron star formations ( jHeger et al.ll2005l ; lOtt et al.ll2006bl ). Besides 
the rapid rotation of the core, convective matter motions and anisotropic neutrino emis- 
sion in the much later postbounce phase are expected to be the primary G W sources with 
comparable amplitudes to the bounce signals. Thus far, various physical ingredients for pro- 
ducing asphericities and the resulting GWs in the p ostbounce phase have been studied, such 
as the roles of pre-collapse density inhomogeneities (IBurrows &: Hayeslll996l ; iMiiler fc Janka 



1 http: / / www.ego-gw.it / 
2 http://geo600. aei.mpg.de/ 



-3 - 



19971 ; lFryerll2004al Jbl) . mod erate rotation of the iron core ( Miiller et al.ll2004j) . nonaxisymmet 



ric rotational instabilities (iRampp et al. 



1998 



Ottetal 



2007b|), g-modes f lOtt et alJl2006af > 
and r-modes pulsations (|Andersson et al.l I2010Q o f protoneutron stars (PNSs), and SASI 



fcotake et alTfl2007L 12009 J bh: iMarek et all J2OO9I ): iMurphv etal] (bo09h ) 



Here SASI, becoming very popular in current supernova researches, is a uni- and bipo- 
lar sloshin g of the stalled super n ova shock with puls a tional strong expansi o n and contrac- 



tion (e. g.. iBlondin et all (120031) : IScheck et al.l (|200J); lOhnishi et al.l (120061 ): iFoglizzo et al. 



( 120071 ): llwakami et al.l (120081 ) ). Based on two-dimensional (2D) simulations that parametrize 
the neutrino heating and cooling by a light-bulb scheme to obtain explosions (which we 
shortly call as parametric SASI simulations) , we pointed out that the GW amplitudes from 
anisotropic neutrino emission increase almost monotonically with time, and that such sig- 
nals m ay be visible to next-generation detectors for a Galactic source (IKotake et al.l 120071 . 
2009 aj). By pe r formin g such a parametric simulation but without the excision inside the PNS, 
Murphy et al.l (120091 ) showed that the GW signals from matter motions can be a good in- 
dicator to get the informati on of the explos i on ge ometry. These features qualitatively agree 
with the ones obtained by lYakunin et al.l (120101 ) who reported exploding 2D simulations 
in which a multi-g roup flux limited diffusion transport is solved with the hydrodynamics. 
Marek et al.l (120091 ) analyzed the GW emission based on their lon g-term 2D (ray-by-ray ) 
Boltzmann simulations, which seem very close to produce explosions (IMarek &: Jankall2009l ). 
They also confirmed that the GWs from neutri nos with cont i nuously growing amplitudes 



(but w ith the different sign of the amplitudes in IKotake et al.l (120071 . l2009al ): lYakunin et al. 
(120101 )). are dominant over the ones from matter motions. They proposed that the third- 



generation class detectors such as the Einstein Telescope are required for detecting the GW 
signals with a good signal-to-noise ratio. 

Thanks to a growing computational power, three-dimensional (3D) simulations, though 
mostly limited to the parametric simulations at present, are now becoming practicable. In 
3D, th e modes of SASI are divided into sloshing modes and spiral modes (e.g.. llwakami et al. 
(120081 )). Asymmetric m = modes so far studied in 2D models and axisymmetric m 7^ 



modes are classified into the sloshing modes, where m stands for the azimuthal index of 
the spherical harmonics Y"/ 71 . In the latter situation, the ±m modes degenerate so that 
the +m modes has the same amplitudes as the — m modes. If random perturbations or 
uniformly rotating flow are imposed on these axisymmetric fl ows in the postbounce phase , 
the degeneracy is b roken and the rotational modes emerge (IBlondin fc Mezzacappal 120071 : 
Iwakami et al.ll2009l ). In this situation, the +m modes has the different amplitudes from — m 
modes. Such rotating non-axisymmetric rn^O modes are called as spiral modes. These non- 
axisymmetric modes are expected to bring about a breakthrough in our supernova theory, 
because they can help to produce explosions more easily compared to 2D due to its extra 
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degree of freedom flNordhaus et al.l 120101). and also because t hey may have a potential t o 
generate pulsar spins (jBlondin fc Mezzacappal (120071 ) , see also lYamasaki &: Fogfizzol ( 120081 ) ; 



Wongwathanarat et al.l (120101 )). 



In a series of our previous papers, we have studied the effects of the sloshing SASI 



modes on th e GW signatures based on the 2D parametric SASI simulations (IKotake et al. 



20071 . l2009al ). In the latter study, we proposed a ray-tracing method to accurately estimate 
GWs generated by anisotropic neutrino emission. Then moving on to the 3D non-rotating 
parametric simulations, we pointed out that the gravitational waveforms in 3D vary much 
more stochastically than for the corresponding 2D models because the explosion anisotropies 
depend sensitively on the growth of the SASI which develops chaotically in all directions 
(IKotake et al.ll2009bl ). As a sequel of these studies, we hope to study the role of the spiral 
SASI modes on the GW emission in this work. We include the effects of rotation as in 



Iwakami et al.l ( 120091 ). That is, when the SASI gradually transits from the linear to nonlinear 
phase, we give a uniform rotation to the flow advecting from the outer boundary of the iron 
core, wh ose specific angul ar momentum is assumed to agree with recent stellar evolution 
models ( jHeger et al.l 120051 ) . After the initial rotational flows advect to the PNS surface, 
rotation begins to have influence over the GW signals due to the non-axisymmetric flow 
motions and anisotropic neutrino emission outside the PNS. Such an approach (excision 
inside the PN S), could be justif ied by some striking evidences that support the slow rotation 
of iron cores (jHeger et al.ll2005l ). and al so that the postbou nce density structure of the PNS 



for such a progenitor (model m!5b6 in iHeger et al 



the ones in the non- rotating models (jOtt et al.ll2006bl ). To see clearly the effects of rotation 



20051 )) has a very similar structure to 



we compute fifteen 3D models in which the initial angular momentum as well as the input 
neutrino luminosities from the PNS are changed in a systematic manner. To estimate the 
GWs from anisotropic neutrino radiation, we perform the ray-tracing analysis. With these 
computations, we hope to clarify how rotation that gives a special direction to the system 
(i.e., the spin axis) , could affect the stochastic GW features obtained in our previous study 



e.g., 



Kotake et al.l (I2009H )) 



The plan of this paper is as follows. In section |2l we shortly summarize the information 
how to construct our 3D models and also how to compute the gravitational waveforms. The 
main results are shown in Sectional We summarize our results and discuss their implications 
in Section HI 
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Numerical Methods and Models 



2.1. Extraction of Gravitational Waveforms 

The numerical met h ods to extract the gravitational waveforms have been already de- 
scribed in iKotake et al.l (l2009al ). We e xtract the GW amp l itudes from mass m otions using 



the standa rd quadrupole formul a (e.g. , IZhuge et al.l (119941 ); iRampp et al.l (Il998f )). and their 



spectra by Flanagan fc Hughesl (119981) with the FFT techniques. F or the GWs generated 
by anisotropic neutrino emission ( IEpsteinlll978l ; iMiiler &: Jankalll997| ). the two modes of the 
GWs can be derived as follows, 



h. 



C [ [ dQ'(l + S (0')c(0'MO + c(6')c(0) x 

JO J4n 

(s(6>)c(<p>)c(0 - c(6')s(0) 2 - s\6')s\<P') dl u (Q',t') 



[s(9')c(4>')c(0 - c(9')s(0} 2 + s 2 (0' W) 



dVL' 



and 



h x = 2C / <m'(l + s(6')c{<J)')s(0 + c(6')c(0) x 

JO Ja-k 

g (^) S (0Q( S (^)c(0Qc(O - c(9')s(Q) dl u (n',t') 
[s{6')c(<f>')c{t) - c(9')s(£)} 2 + s 2 (6>)s 2 (<j)') dW ' U 

where s(A) = sin(A),c(B) = cosB, C = 2G/(c 4 R) with G, c and R, being the gravitational 
constant, the speed of light, the distance of the source to the observer respectively, dl u /dQ 
represents the direction-dependent neutrino lum inosity emitte d per u nit of solid angle into 
direction of Q, and £ is the viewing angle (e.g., IKotake et al.l (l2009al )). For simplicity, we 
consider here two cases, in which the observer is assumed to be situated along 'polar' (£ = 0) 
or 'equatorial' (£ = vr/2) direction. It is important to note that in the case of 2D axisymmetric 
case, the only non-vanishing component is the plus mode for the equatorial observer, 

dl u (6',t') 



hi 



2C I di 
'o 



60 



dVL' 



(3) 



where the function of $(#') nas positive values in the north polar cap for < 9' < 60° 
and in the south polar cap for 120° < 6' < 180°, but becomes negative values between 
60° < 6' < 120° (see Figure 1 of IKotake etHl f j2007h ). 



To determine dl v /dVt in equations dU [2]), we perform a ray- tracing calculation which 
we shortly summarize in the following (see IKotake et al.l (j2009al ) for more detail). In the 
ray-tracing approach, we consider transfer along the ray specified by a constant impact 
parameter p. The coordinate along p is called s, satisfying 

r = (p 2 + s 2 ) l '\ (4) 
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where r is the radial coordinate. In order to get the numerical convergence of dl v (fl)/dQ, 
we need to set 45,000 rays for each direction, which consists of 500 x 90 rays, where the 
former is for the impact parameters covering from the inner- (p in = 50 km) to the outer- 
boundary (pout — 2000 km) of the computational domain and the latter is for covering the 
circumference (e.g., 2tt) of the concentric circles on the plane perpendicular to the rays. 

The transfer equation of the neutrino occupation probability f u (e u ,p,s) for a given 
neutrino energy e v along each ray is given by, 

j- = j{e v ,p,s)(l - fv(e v ,p,s)) , (5) 

where j and A is the emissivity and abs o rptivi t y via neutrino absorp tions and emission by 



free nucleons {y e + n e + p) ( iBruennl ( 119851 ): lOhnishi et al.l (120061 ) ) . which are dominant 



processes outside the PNSs. The optical depth for those reactions are estimated by r v = 
J r °° 1/A. For the sake of simplicity, the neutrino scattering and the velocity-dependent terms 
in the transport equation are neglected. Along each ray, f v is transferred by the line integral. 
When the line integral starts from the surface on the PNS, we set the initial value of 

f( e ") = i - T // rp \ ■ j-> (6) 

1 + exp{e u /k B T u ) 4tt 

assuming that the neutrino distribution function at the surface is approximated by the Fermi- 
Dirac distribution with a vanishing chemical potential. Here the neutrino temperature is set 
to be constant near T Uc = 4 MeV, whose values change slightly depending on the input 
neutrino luminosity. Note that these values are constant in ti me for each m o del. This is 



necessary to realize the steady unperturbed initial states (e.g., lOhnishi et al.l (120061 ) ) . For 
the rays that do not hit the PNS, we start the line integral from the outer most boundary 
antipodal to the line of sight, where /„ is essentially zero. 

By a post process, we perform the line integral up to the outer-most boundary for 
each hydro-timestep. With f(e u ,p,s out ), which is obtained by the line integral up to the 
outer-most boundary, the neutrino energy fluxes along a specified direction of ft can be 
estimated, 

d '" in - p) ' /(^o„.) • («J • (7) 



dads J ( 2l fc)3- 

By summing up the energy fluxes with the weight of the area in the plane perpendicular to 
the rays, we can find dl u /dVL along a specified direction O, 



dVt J dVtdS L. dfldS 



pi, 



Repeating the above procedures, dl u (Q) / dQ can be estimated for all the directions, by which 
we can find the amplitudes of the neutrino GWs through equations (fTEZl) . In the following 
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computations, we assume that the distance to the GW source is comparable to our galactic 
center (R = 10 kpc) unless stated otherwise. 



2.2. Construction of 3D Models with Rotation 



The employed nu merical methods and model concepts are essentially the same as those 
in our previous paper ( jlwakami et al.l 120091 ). Using the ZEUS-MP code (jHayes et al.ll2006l ) 
as a hydro-solver, we solve the dynamics of the standing accretion shock flows of matter at- 
tracted by the protoneutron star and irradiated by ne utrinos emitted f r om th e PNS. We em- 



ploy t he so-called light-bulb approximation (see, e.g.. Ijanka fc Mullerl (119961 ); lOhnishi et al. 



(120061 )) and adjust the neutrino luminosities from the PNSs to trigger explosions. This 
method makes it possible for us to study the properties of GWs in the postbounce phase, 
namely from th e shoc k- stall, through the grow t h of the SASI, to 3D explosions (see also 
Nordhaus et al.l ( 120101 ): I Wongwathanarat et al.l ( 120101 )). and is found to work well in 2D, 
capturing th e essential f eatures obtained by more realistic simulations (see references in 
Scheck et al.l tool S3)). 



The computational grid is comprised of 300 logarithmically spaced, radial zones to cover 
from the absorbing inner boundary of ~ 50 km to the outer boundary of 2000 km, and 30 
polar (6) and 60 azimuthal (<p) uniform m esh points, wh i ch are used to cover the whole 
solid angle (see for the resolution t ests in llwakami et al.l ( 120081 )). The initial conditions 
are provided in the same manner of lOhnishi et al.l tood \. whi ch describes the sphericall y 
symmetric steady accretion flow through a standing shock wave (lYamasaki &: Foglizzoll2008l ). 



In constructing the initial conditions, we assume a fixed density p- m = 10 11 g cm -3 at the 
inner boundary. And the initial mass accretion rates and the initial mass of the central 
object are set to be M = 1 M s -1 and M in = 1.4 M , respectively. To include rotation, we 
impose a rigid rotation on the outer boundary of the computational domain as follows, 



= ^vl D (r) sm6, (9) 

where v^ D denotes the unperturbed <p component of velocity, v l r D (r) — —M /47rr 2 p(r) is the 
unperturbed radial velocity, and (3^ denotes the rotation parameter (see llwakami et al.l (120091 ) 
for more detail). We examine the flow characteristics in the following three ways = 0, 0.01, 
and 0.015 which corresponds to the specific angular momentum L ~ (4 — 6) x 10 15 cm 2 s~ x 
on the equatorial plane. This choice is close to the recent rotating presupernoya mo dels 
that include the effects of magnetic breaking during stellar evolution ( jHeger et al.ll2005l ). It 
is computationally prohibitive at present to follow a long-term postbounce evolution in 3D 
simulations that implement the spherical coordinates due to its severe Courant-Friedlichs-Lax 
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condition around the pole. As a prelude to the full 3D simulations starting from gravitational 
collapse to explosions without excision inside the PNS, we choose to take into account the 
effect of rotation by the exploratory method as mentioned above. To induce non-spherical 
instability we add random velocity perturbations to be less than 1 % of the unperturbed 
radial velocity. By changing also the neutrino temperature at the PNS surface (e.g., equation 
(6)), we compute fifteen 3D models whose input neutrino luminosities vary in the range of 
L Ve = 6.0 - 6.8 x 10 52 erg s" 1 (see Tabled]). 

In Tabled], At represents the simulation time when the average shock radius that contin- 
uously increases with the growth of SASI, reaches the outer boundary of the computational 
domain with a typical explosion energy of ~ 10 51 erg. On the other hand, for models denoted 
by " — " , we terminated the simulation at about 1000 ms, not seeing the increase of the shock 
radius (represented by "No" in the table). It can be seen that models with higher neutrino 
luminosities such as models of series D and E produce explosions without rotation (models 
DO and E0), while models with the intermediate luminosities such as models B2, CI and C2 
produce explosions only when they possess rotation (compare models BO and CO). This is 
mainly because the centrifugal force induced by rotation makes not only the shock radius 
but also the gain region larger, which works to assist the onset of explosion. Regardless of 
rotation, no explosions are obtained in the lowest luminosity models (model of series A) for 
the rotation parameters investigated here. 

To analyze clearly the role of rotation on the GW emission, models of series A that 
do not produce explosions, are most favorable. As will be explained more in detail later, 
this is because in the higher luminosity models, the violent fluid motions as well as vigorous 
convective overturns, which works to smear out the effects of rotation, make the analysis 
more complicated. For example, At becomes shorter with larger angular momentum for the 
intermediate luminosity (models of series B and C in Tabled]). However this trend is no more 
true for models of series D (e.g., model D2). This non-monotonic effect can be understood 
as a consequence of stochasticity in a regime of high neutrino luminosity. 

a s k in Table 1 represents the angle between spin and kick direction s that the central PNS 



may r eceive at the final simulation time, which we estimate according to lWongwathanarat et al 



(120101 ). The obtained spin-kick angle ranges from 25° (model B2) in which the spin direction 
is marginally aligned with the kick direction, 88° (model E0) in which the two directions 
are closely perpendicular, to 170° (model CI) in which the two directions are closely op- 
posite. When a rapid rotation is imposed in our models (e.g., for our models of series 2), 
the final spin direction is closely aligned with the direction of the induced rotation (namely, 
perpendicular to the equatorial plane), while even in this case, stochasticity of the kick di- 
rection seems rather unaffected. To clarify the correlation clearly, a more elaborate study 
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is needed in which the initial neutrino luminosity, perturbations, and rotation, are varied 
much more systematically, which we plan to study as a sequel of this study (Iwakami et al. 
in preparation). 

As already mentioned, we can adjust the epoch by ha nd when the rotational flow advects 



to the PNS surface (see Figure 1 in llwakami et al.l (120091 ) ). The deformation of the standing 



shock becomes remarkable typically at about t = 100 ms, marking the epoch when the SASI 
transits from the linear to non-linear regime. For the non-exploding models of Al and A2, 
we choose to inject the rotational flow so that it advects to the PNS surface at around 
t = 400 ms when the SASI has been long developed into the non-linear phase. For the other 
luminosity models, we take the time at around t = 100 ms (therefore soon after the SASI 
enters the non-linear phase). This is because some models with higher neutrino luminosities 
show a trend of explosions before the rotational flows advect to the central regions if we 
delay the epoch of the injection and also because we want to fix the epoch common to these 
models. 
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Table 1. Model Summary 



Model 




A* 


At 


Explosion? 


l/i po1 

pmax | 


\ "max | 


l/i po1 1 

| ,t i^,max | 


^i/,max | 








(10 52 erg/s) 




(ms) 




(10- 22 ) 


(10-22) 


(lO- 22 ) 


(lO" 22 ) 


(IO-^Mqc 2 ) 


(°) 


AO 


6.0 







No 


0.98(+) 


1.59(x) 


0.72(+) 


1.13(x) 


2.54 


100 


Al 


6.0 


0.01 




No 


1.98(+) 


7.62(+) 


1.18(+) 


6.74(+) 


6.76 


104 


A2 


6.0 


0.015 




No 


2.58(+) 


16.7(+) 


0.91(x) 


15.1(+) 


13.4 


63 


BO 


6.2 







No 


2.93(+) 


3.19(+) 


2.45(+) 


2.38(+) 


5.65 


68 


Bl 


6.2 


0.01 




No 


2.88(x) 


11.7(+) 


1.37(x) 


10.3(+) 


15.1 


55 


B2 


6.2 


0.015 


721 


Yes 


5.03(x) 


11.2(+) 


2.86(x) 


5.90(+) 


16.1 


25 


CO 


6.4 







No 


4.11(+) 


5.00(+) 


3.46(+) 


4.28(+) 


10.7 


41 


CI 


6.4 


0.01 


650 


Yes 


4.62(+) 


6.46(+) 


2.05(x) 


4.97(+) 


12.1 


170 


C2 


6.4 


0.015 


515 


Yes 


3.43(+) 


9.38(+) 


l.ll(x) 


6.09(+) 


8.96 


141 


DO 


6.6 





550 


Yes 


4.31(+) 


2.93(+) 


2.25(+) 


2.14(x) 


8.18 


113 


Dl 


6.6 


0.01 


460 


Yes 


4.73(x) 


4.70(x) 


3.72(x) 


3.91(x) 


6.17 


45 


D2 


6.6 


0.015 


525 


Yes 


2.31(x) 


4.20(+) 


1.07(x) 


3.50(+) 


6.87 


82 


E0 


6.8 





510 


Yes 


3.44(x) 


2.56(x) 


2.03(x) 


2.06(x) 


6.60 


88 


El 


6.8 


0.01 


425 


Yes 


2.94(x) 


2.33(x) 


0.94(+) 


0.96(x) 


6.67 


150 


E2 


6.8 


0.015 


410 


Yes 


2.34(x) 


3.49(x) 


9.69(x) 


1.69(+) 


5.44 


168 



Note. — L„ c denotes the input luminosity. At represents the simulation time. 0$ is a rotational parameter introduced in 
equation (9). /tm'ax, C represents the maximum amplitudes (neutrino + matter) during the simulation time seen from the 
pole or the equator, respectively, while /i p ^ ax , /i^ ax are the ones including only the neutrino contribution. The polarization 
of the GWs at the maximum is indicated by + and X. Eqw, is the radiated energy in the form of the neutrino GWs in unit 
of Mqc 2 . a s k is the spin-kick angle (see text for more detail). The distance to the source is assumed to be 10 kpc. 
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3. Results 

First of all, we summarize how the inclusion of rotation could have impacts on the 
hydrodynamics as well as the overall trends in the GW emission in section 13.11 Then in 
section I3.2[ we move on to analyze the gravitational waveforms more in detail. 



3.1. Stochastic Hydrodynamics and GWs 

Figures [Tf2l show the 3D hydrodynamic features of the SASI near at the shock breakout 
from the outer boundary of the computational domain. Comparing the top left (model B2) 
to the bottom left panel (model C2) in Figure [T], it can be inferred that the direction of 
explosion varies models to models. Given the same initial angular momentum, the major 
axis of the growth of SASI happens to be toward the rotational axis for model B2 (top left 
panel showing the propagation of the high entropy bubbles (colored by red) closely parallel 
to the rotational axis (:2-axis)), while it is shown to be toward the equatorial direction (:y- 
axis) for model C2 (bottom left). A common feature seen for the pair models is the spiral 
flows inside the shock that rotates globally around the z-axis, which is seen like arcs (right 
panels). This is as a result of introducing the rotation whose axis corresponded to the z-axis. 
As a faster rotation was added to the flows, spiral flows run with higher velocity from the 
triple points that are formed inside the standing shock, to near the PNS, and then flows 
rotate with higher velocity inside the high-entropy blobs. 

Figure [2] shows the hydrodynamic features for models with higher neutrino luminosities 
(models DO (top) and E2 (bottom)). Due to the higher neutrino luminosity, the volumes 
of high entropy blobs (colored by red) behind the shock become much larger than the ones 
in Figured] (right panels). Some spiral structures are depicted in the bottom right panel of 
Figure [21 however, it is more difficult to see the structures compared to Figure 1. Observing 
from the horizontal direction, model DO explodes in a rather spherical manner (top left), 
while the major axis of SASI for model E2 is tilted about 45 degrees from the rotational axis 
(bottom left). Our 3D results presented here imply that the direction of explosion as well 
as the major axis of the growth of SASI cannot be predicted a priori, reflecting the chaotic 
fluid motions driven by the interplay of the neutrino heating, SASI, and co nvection. This 



feature is qualitatively sam e as the previous 3D results bu t without rotation (jlwakami et al. 



20081 : iNordhaus et al.ll2010l : IWonewathanarat et al.ll2010h . 



As the explosion dynamics becomes chaotic, there is no systematic dependence of the 
maximum GW amplitudes and the radiated GW energies on the input neutrino luminosities 
(see Table [1]). In fact, the largest emitted GW energy is obtained for model B2 with the 
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intermediate neutrino luminosity (see, i?GW Tabled]). Regarding the polarization (+ or x) 
of the maximum GW amplitudes, it may be hard to decipher any characteristics in Table [I] 
due to its messiness. But looking very carefully, one may be able to find that the polarization 
of the maximum GWs from neutrinos seen from the equator (/i^ax) * s always plus (+). 

This feature can be more clearly seen in the waveform catalogues presented from Figure 
[3] to [6l From Figure 3, it can be seen that h e ^_ (blue line) becomes prominent among others 
after around t ~ 500 — 600 ms for models with rotation (middle and bottom panels). This 
trend seems quite general as seen in the other luminosity models (see middle and bottom 
panels in Figure S] to [6]). Note again that the time when the rotational flow is adjusted 
to touch the PNS surface is t = 400 ms for models of series A (and t = 100 ms for the 
other) , which accounts for the time difference in the rise of h e ^ (compare Figure 3 to Figure 
4 to 7). In the absence of rotation, the waveforms do not exhibit such a feature and vary 
much more stochastically (e.g., top panels in Figure [3] to [7]) . Comparing the bottom panels 
in Figures [3] to [7J it can be seen that the increasing trend becomes more fainter as the 
input luminosity becomes larger. This is mainly because vigorous matter motions irradiated 
by intense neutrino heating smear out the growth of the spiral SASI modes as well as the 
structures of the spiral flows. To extract purely the rotational effects on the GWs, it is most 
convenient to analyze models of series A that fail to produce explosions. In the following, we 
take models of series A as a reference and look more in detail the reason why the stochastic 
nature of the GWs is broken (or weakened) due to the inclusion of rotation. 
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Fig. 1. — Snapshots of the entropy distributions for models B2 (top panels, t = 624 ms) and 
C2 (bottom panels, t = 513 ms) seen from the equator (left panels, = 0) or the pole (right 
panels, 9 = 0), respectively. The second and fourth quadrant of each panel shows the surface 
of the standing shock wave. In the first and third quadrant, the profiles of the high entropy 
bubbles (colored by red) inside the section cut by the ZY (left) or YX (right) plane are 
shown. The side length of each plot is 1000km. The insets show the gravitational waveforms 
with '+' on each curves representing the time of the snapshot. Note that the colors of the 
curves are taken to be the same as the top panel of Figure [31 
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Fig. 3. — Gravitational waveforms from the sum of neutrinos and matter motions (left) and 
only from neutrinos (right) for models AO (top), Al (middle), and A2 (bottom). The time 
is measured from the epoch when the neutrino luminosity is injected from the surface of 
the neutrino sphere. In all the computed models, SASI gradually transits from the linear 
to non-linear regime at about 100 ms, simultaneously making the amplitudes deviate from 
zero. For models of series A shown here, the rotational flow is adjusted to advect to the 
PNS surface at around t = 400 ms (see text for more detail). The supernova is assumed to 
be located at the distance of 10 kpc. 
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Fig. 4. — Same as Figure [3] but for models of series B . Note that except for models of series 
A, the rotational flows are adjusted to advect to the PNS surface at around t = 100 ms. 
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Fig. 5. — Same as Figure H] but for models of series C. 
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Same as Figure H] but for models of series D. 



- 19 - 




Pole, + 
Pole.x 
Equator,+ 
-3 Equator.x 







100 



200 300 
Time [ms] 



400 



500 



4 

3 

„ 2 

OJ 
CM 

a 1 

O 

E 

<D 
c 

t -1 

-2 
-3 



Model E0 (L v =6.8x1 52 erg/s, no rotation) 




Pole, + 
Pole.x 
Equator,+ 
Equator.x 



100 



200 300 
Time [ms] 



400 



500 



t -1 
-2 
-3 



Model E1 (L v =6.8x1 52 erg/s, moderate rotation) 




W I' 



nil 



Pole, + 
Pole.x 
Equator,+ 
Equator.x 



50 100 150 200 250 300 350 400 
Time [ms] 



4 - 
3 - 
2 - 
1 - 


-1 - 
-2 



Model E2 (L v =6.8x1 52 erg/s, rapid rotation) 



4 



Pole, + 
Pole.x 
3 Equator,+ 
Equator, x 




-4 



50 100 150 200 250 300 350 400 
Time [ms] 



o 

^ 



Model E1 (Ly =6.8x1 52 erg/s, moderate rotation) 




Pole, + 
Pole,x 
Equator,+ 
Equator,x 



50 100 150 200 250 300 350 400 
Time [ms] 



Model E2 (L v =6.8x1 52 erg/s, rapid rotation) 




Pole, + 
Pole.x 
Equator,+ 
Equator.x 



50 100 150 200 250 300 350 400 
Time [ms] 



Fig. 7. — Same as Figure H] but for models of series E. 
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3.2. Breaking of Stochasticity Due to Rotation 

Figure [8] illustrates two typical snapshots of the flow fields for model A2 before and 
after the rotational flow approaches to the PNS surface (left and right panels), seen from 
the pole (top) or from the equator (bottom), respectively. The left panels indicate that the 
deformation of the shock surface is only mild for this model, although the SASI has already 
entered the non-linear phase at around t = 100 ms. From the bottom right panel, one may 
guess the presence of the sloshing modes that happen to develop along the rotational axis 
(z-axis) at this epoch. It should be emphasized that the dominance of observed in the 
current 3D simulatio ns have nothing to do with the one found in our previous 2D studies 



( jKotake et al.ll2009al ). Free from the 2D axis effects, the major axis of the SASI changes 
stochastically with time, and the flow patters behind the standing shock simultaneously 
change in every direction. As a result, the sloshing modes can make only a small contribution 
to the GW emission (e.g., the inset in the left panels of Figure [8]). The remaining possibility 
is that the spiral flows clearly seen in the top right panel should be a key importance to 
understand the GW feature that we pointed out in the last section. 

The left panel of Figure [9] shows the local neutrino energy fluxes {:dl u / (dQdS), equation 
(J7J)) seen from the northern hemisphere. Seen as bright arcs in the left panel, the higher 
values of the neutrino energy fluxes are shown to coincide with the high temperature regions 
in the right panel. This is because the compression of matter is more enhanced in the vicinity 
of the equatorial plane due to the presence of the spiral flows (see velocity vectors in the 
right panel). 

Figure [TU] shows the time evolution of dl u /dQ (equation OH])) in the vicinity of the north 
pole (9 = 0), the equator (9 = 7r/2), and the south pole (9 = ir) (left panels), and their 
differences from the equator (right panels) for models A2 (top panels) and AO (bottom 
panels), respectively. A common feature seen from the right panels is that the dominance 
of the neutrino emission in the north (red line) and south poles (green line) is occasionally 
anti-correlated. This is a consequence of the low-mode nature of SASI, here of £ = 1. The 
correlation shown here as well a s the resulting GWs is much weaker than the ones in our 



2D study ( IKotake et al.l ( j2009a[ )). because the major axis of the SASI changes much more 
randomly. 

The most important message in Figure [TU] is the gradual deviation of dl v j dfl for model 
A2 seen in the top right panel. This feature becomes remarkable only after t ~ 500 — 600 
ms, because it takes ~ 100 — 200 ms for the spiral SASI modes to develop behind the 
standing shock after the spiral flows approached to the PNS surface (:t = 400 ms denoted 
by the vertical line). As already mentioned, the dominance of the neutrino luminosity seen 
from the polar directions is due to the spiral flows that develop near in the vicinity of 
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the equatorial plane. Comparing the two panels in Figure HU a contrast of the neutrino 
luminosity in the polar regions (9 ~ 0° or 9 ~ 180°) to other regions becomes prominent 
only after t ~ 600 ms for model A2 (right panel), which is visualized as stripes bridging 
between the north and south pole. Without rotation, the anisotropy of neutrino radiation 
remains to be much smaller (left panel in Figure ITT]) , leading to much smaller GW amplitudes. 
Remembering again that in equation fl3]) is positive near the north and south polar 

caps, the dominance of the polar neutrino luminosities leads to make the positively growing 
feature in h e ^+, so far depicted in models that include rotation (e.g., Figures [3}17j) . It is noted 
here that the lateral-angle (9) dependence of equation (1) can be approximated by equation 
(3) because the neutrino anisotropy in the azimuthal direction is much weaker than the one 
in the lateral direction. 

Figure d2] shows how the above GW features could or could not change if we change the 
strength of the initial velocity perturbation (top panels), the epoch when the initial rotational 
flow advects to the center (bottom left panel), and the numerical resolution (bottom right 
panel), respectively. In the top two panels, larger initial (velocity-) perturbations (5%) are 
imposed (for model Al (left) and A2 (right) as indicated by pert5%) in contrast to the 
fiducial value of 1% (e.g., section 2.2). Due to the large perturbation, the timescale when 
the non-linear phase sets in, becomes typically 30 — 40 ms earlier compared to the fiducial 
models (compare the middle and bottom panels of Figure 3 to the top two panels of Figure 
[j~2l . Regardless of the difference, the increasing trend of the GWs (blue lines) is shown to 
be unaffected. This is also the case when the rotational flow is adjusted to advect to the 
PNS in the linear SASI phase (bottom left panel), and when a finer numerical resolution 
is taken (bottom right panel). In the bottom right panel, the numerical resolution for the 
azimuthal direction is doubled compared to the fiducial value of 60 mesh points (indicated 
by "high-res" for model A2). For the bottom left panel, the initial perturbation is taken to 
be as small as 0.1% to make the linear phase longer. For this model, the rotational flow 
advects to the PNS about t ~ 100 ms in the linear SASI phase, and the transition to the 
non-linear phase takes place at t > 200 ms. As is shown, the increasing trend appears also 
from the linear SASI phase (t ~ 100 ms). 

Finally Figure [T31 depicts the GW spectra for model AO (left panel) and A2 (right panel). 
The neutrino GWs seen from the equator (green line, right panel) for model A2 is larger 
than the one seen from the pole, and it is slightly larger compared to model AO in the 
lower frequencies below ~ 10 Hz. It is true that the GW signals from neutrinos are very 
difficult to detect for ground-based d etectors whose sensitivity is limited mainly by the seis 



mic noises at such lower frequencies (lAndo fc the TAMA Collaboration! 120051 ; lAbbott et al. 



2005 ; Weinstein 2002 ; Kuroda fc LCGT Collaboration 201ol ). However these signals may be 



detectable by the recently proposed future space interferometers like Fabry-Perot type DE- 
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CIGO (|Kawamura et al.l ( 120061 ). black line in Figure IT3|) . The sensitivity curve is taken from 
Kudoh et al.l ( 120061 ). Contributed by the neutrino GWs in the lower frequency domains, the 
total GW spectrum tends to become rather flat over a broad frequency range below ~ 100 
Hz. These GW features obtained in the context of the SASI- aided neutrino-driven mech- 
anism are different from the ones expected in the other cand i date supernova mechanism , 
such as the MHD mechanism ( e.g.. lObergaulinger et al.l (120061 ) ; iTakiwaki fc Kotakd (120101 ) ) 



and the acoustic mechanism ( jOtt et al.ll2006al ). Therefore the detection of such signals is 



expected to provide an important probe into the long-veiled explosion mechanism. 

Regarding the GW amplitudes from matter motions, we find no significant change or 
increase due to the spiral flows either seen from the polar or equatorial directions. This is 
mainly because the spiral flows can develop only outside the PNS taken to be 50 km in our 
idealized simulations, leading to have small changes in the mass-quadrupole (due to its small 
masses) . Thi s situa tion might be akin to the one observed in the early 3D simulations by 
Rampp et al.l ( 119981 ). To produce sizable GW amplitudes, the non-axisymmetric i nstabilities 



that globally develop in the more central regions to the PNS should be needed (e.g.. lOtt et al. 



( 12007bl ); IScheidegger et al.l ( 12010I )). To explore these fascinating phenomena, one apparently 
needs to perform full 3D radiation-hydrodynamic simulations covering the entire stellar core 
and starting from gravitational collapse to explosions in a consistent manner. 
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Fig. 8. — Partial cutaway of the entropy isosurfaces and the velocity vectors on the cutting 
plane for model A2. Left and right panels are at t = 300 ms and t = 900 ms corresponding 
to the epoch before and after the rotational flow approaches to the PNS surface, respec- 
tively. Top and bottom panels are for the polar and equatorial observer, respectively. The 
insets show the gravitational waveforms with '+' on each curves representing the time of the 
snapshot. Note that the colors of the curves are taken to be the same as the top panel of 
Figure El 
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Fig. 9. — Left and right panels show the neutrino energy fluxes of dl u /(dfldS) (equation(J7|)) 
seen from the northern hemisphere and the temperature distributions in the equatorial plane 
for model A2 at t = 900 ms. 




Time [s] Time [s] 

Fig. 10. — Time evolution of the directional dependent neutrino luminosity: dl u /dQ (left 
panels) in the vicinity of the north pole, the equator, and the south pole, and their differences 
from the equator (right panels) for models A2 (top panels) and AO (bottom panels). Vertical 
lines in the right panels represent the epoch of t = 400 ms when the rotational flow advects 
to the PNS surface. 
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Fig. 11. — Contour of dl v /dVL ((azimuthal) angle-averaged) as a function of time (horizontal 
axis) and the polar angle (vertical axis) for models AO (left panel) and A2 (right panel). 
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Fig. 12. — Variations of model series of Al and A2. For top panels, larger initial (velocity-) 
perturbations (5%) are imposed (for model Al (left) and A2 (right) indicated by pert5%) in 
contrast to the fiducial value of 1% (e.g., section 2.2). For the bottom left panel, the initial 
rotational flow is tuned to advect to the PNS surface in the linear SASI phase for model 
A2 (indicated by A2-linear). In the bottom right panel, the numerical resolution for the 
azimuthal direction is doubled compared to the fiducial value of 60 mesh points (indicated 
by "high- res" for model A2). 
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Fig. 13. — Spectral distributions of GWs from matter motions ("Matter") and neutrino 
emission ( "Neutrino" ) seen from the pole or the equator for models AO (left panel) and A2 
(righ t panel) with the expected detecti on limits of TAMA dAndo fc the TAMA Collaboration 
2005h . first LIGO JAbbott et al.ll2005h . advance d LIGO JWeinsteinll2002h. Large-scale Cry~ 
genie Gravitational wave Te lescope (LCGT) (IKuroda fc LCGT Collaboration! 120101 ) and 



Fabry-Perot type DECIGO feawamura et all l200fil : iKudoh et all |20oJ The distance to 
the supernova is assumed to be 10 kpc. Note that for the matter signal, the + mode seen 
from the polar direction is plotted. 
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4. Summary and Discussion 



We studied how the spiral modes of the SASI can have impacts on the properties of GWs 
by performing 3D simulations that mimic the SASI-aided core-collapse supernovae. To see 
the effects of rotation, we impos ed a uniform ro t ation on the flow advecting from the outer 
boundary of the iron core (as in llwakami et al.l (120091 )). whose specific angular momentum 
is assumed to agree with recent stellar evolution models. We computed fifteen 3D models in 
which the initial angular momentum as well as the input neutrino luminosities from the PNS 
are changed in a systematic manner. By performing a ray-tracing analysis, we accurately 
estimated the GW amplitudes from anisotropic neutrino emission. With these computations, 
we found that the gravitational waveforms from neutrinos in models that include rotation 
exhibit a common feature otherwise they vary much more stochastically in the absence of 
rotation. The breaking of the stochasticity stems from the excess of the neutrino emission 
parallel to the spin axis. This is because the compression of matter is more enhanced in 
the vicinity of the equatorial plane due to the growth of the spiral SASI modes, leading to 
the formation of the spiral flows with higher temperatures circulating around the spin axis. 
We pointed out that a recently proposed future space interferometers like Fabry-Perot type 
DECIGO would permit the detection of these signals for a Galactic supernova. 

It should be noted that the approximations taken in the simulation, such as the excision 
inside the PNS with its fixed inner boundary and the light bulb approach with the isentropic 
luminosity constant with time, are only a very first step towards realistic 3D supernova 
simulations. As already mentioned, the excision of the central regions inside PNSs hin- 



ders th e efficient gravitational emission there, s uch as by the g-mode oscillations (lOtt et al. 



2006al ) and the non-axisymmetric instabilities (IQtt et al.ll2007bl ; IScheidegger et a 



the PNSs, and the enhanced neutrino emissions inside the PNSs (IMarek et al 



should also affect the peak frequency of the GW spectra as pointed out by iMurphy et al. 



201of) of 
2009h . It 



(120091 ). The rotational flows advecting into the inner core should spin up the PNS, which 
is expected to conversely affect the dynamics outside (IRantsiou et al.l 120101 ) . These impor- 
tant feedback should be taken into account to unravel the effect of rotation more precisely. 
Remembering these caveats in mind, one encouraging news to us is that the gravitati onal 
waveforms obtained in the 2D radiation-hydrod ynamic simulations (jYakunin et al.ll2010l ) are 
similar to the ones obtained in our 2D study (IKotake et al.l l2009al ). It is also noted that 



magnetic effects sh ould have impacts not only on the growth of the SASI (jEndeve et al.ll2010 



Guilet et al.ll2010l ). but also on the GW waveforms, which are remained to be studied in the 
3D MHD simulations. Nowadays, the next generation 3D results are being reported by using 
new simulation techniques such as th e Yin- Yang grids (IWongwathanarat et al.ll2010l ) and the 
adaptive mesh refinement approach (INordhaus et al.l 120101 ) . We hope that our results may 
give a momentum to theorists for making the GW prediction with a better quantitative 
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precision based on the sophisticated 3D supernova modeling. 
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